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Abstract 

We consider pure-jump transaction-level models for asset prices in continuous 
time, driven by point processes. In a bivariate model that admits cointegration, we 
allow for time deformations to account for such effects as intraday seasonal patterns 
in volatility, and non-trading periods that may be different for the two assets. We 
also allow for asymmetries (leverage effects). We obtain the asymptotic distribution 
of the log-price process. We also obtain the asymptotic distribution of the ordinary 
least-squares estimator of the cointegrating parameter based on data sampled from 
an equally-spaced discretization of calendar time, in the case of weak fractional coin- 
tegration. For this same case, we obtain the asymptotic distribution for a tapered 
estimator under more general assumptions. In the strong fractional cointegration 
case, we obtain the limiting distribution of a continuously-averaged tapered estima- 
tor as well as other estimators of the cointegrating parameter, and find that the rate 
of convergence can be affected by properties of intertrade durations. In particular, 
the persistence of durations (hence of volatility) can affect the degree of cointegration. 
We also obtain the rate of convergence of several estimators of the cointegrating pa- 
rameter in the standard cointegration case. Finally, we consider the properties of the 
ordinary least squares estimator of the regression parameter in a spurious regression, 
i.e., in the absence of cointegration. 
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1 Introduction 



The increasingly widespread availability of transaction-level financial price data motivates 
the development of models to describe such data, as well as theory for widely-used statistics 
of interest under the assumption of a given transaction-level generating mechanism. We 
focus here on a bivariate pure-jump model in continuous time for log prices proposed 
by Hurvich and Wang (2009, 2010) which yields fractional or standard cointegration. The 
motivation for using a pure-jump model is that observed price series are step functions, since 
no change is possible in observed prices during time periods when there are no transactions. 
Examples of data sets that would fit into the framework of this model include: buy prices 
and sell prices of a single stock; prices of two different stocks within the same industry; 
stock and option prices of a given company; option prices on a given stock with different 
degrees of maturity or moneyness; corporate bond prices at different maturities for a given 
company; Treasury bond prices at different maturities. 

Though our paper is not entirely focused on the case of fractional cointegration, we 
present here some evidence that this case may arise in practice in financial econometrics. 
We considered option and underlying best-available bid prices for 69 different options on 
IBM at 390 one-minute intervals from 9:30 AM to 4 PM on May 31, 2007. Using a log- 
periodogram estimator based on 390°'^ frequencies, we found that the logs of the original 
series had estimated memory parameters close to 1, while the residuals from the OLS re- 
gression of the log stock price on the log option price had estimated memory parameters 
that were typically less than 1. Specifically, of the 69 estimated memory parameters based 
on these residuals, the values ranged from 0.05 to 1.14 with a mean of 0.55 and a stan- 
dard deviation of 0.28, with 30 of these estimates lying between 0.5 and 1, while 32 were 
between and 0.5. Thus, there is evidence for cointegration in most of the series studied, 
and often the evidence points towards fractional rather than standard cointegration. Fur- 
thermore, the OLS estimate of the cointegrating parameter (assuming that cointegration 
exists) ranged from —0.21 to 0.39, with a mean of 0.04 and a standard deviation of 0.13. 
This provides evidence that the cointegrating parameter is in general not equal to one in 
the present context, so it is of interest to study properties of estimates of this parameter. 

Two basic questions that we address in this paper are the asymptotic distribution of 
the log prices as time t — >■ oo, and of the usual OLS estimator of the cointegrating param- 
eter based on n observations of the log prices at equally-spaced time intervals as n — )■ oo. 
Most of the existing methods for deriving such limit laws (see Robinson and Marinucci 
(2001)) cannot be applied here because the continuous-time log- price series are not diffu- 
sions and because the discretized log-price series are not linear in either an iid sequence, a 
martingale difference sequence or a strong mixing sequence. Nevertheless, it is of interest 
to know whether and under what conditions the existing limit laws, based, say, on lin- 
earity assumptions in discrete time, continue to hold under a transaction-level generating 
mechanism. 

In the model of Hurvich and Wang (2010, 2009) the price process in continuous time is 
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specified by a counting process giving the cumulative number of transactions up to time t, 
togetlier witli tlie process of changes in log price at the transaction times. This structure 
corresponds to the fact that most transaction-level data consists of a time stamp giving 
the transaction time as well as a price at that time. In such a setting, another observable 
quantity of interest is the durations, i.e., the waiting times between successive transactions 
of a given asset. There is a growing literature on univariate models for durations, including 
the seminal paper of Engle and Russell (1998) on the autoregressive conditional duration 
models (ACD), as well as Bauwens and Veredas (2004) on the stochastic duration model 
(SCD), and Deo ct al. (2010) on the long-memory stochastic duration model (LMSD). 

Deo ct al. (2009) showed that, subject to regularity conditions, if partial sums of cen- 
tered durations, scaled by n"*^"'"'"^/^) with d G [0,1/2), satisfy a functional central limit 
theorem then the counting process N{t) has long or short memory (for d > 0, d = 0, 
respectively) in the sense that VarA^(t) ~ Ct'^'^'^^ as t — !■ oo (with C > 0), and they gave 
conditions under which this scaling would lead to long memory in volatility. In particular, 
LMSD durations with d > lead to long memory in volatility. The latter property has 
been widely observed in the econometrics literature, while evidence for long memory in 
durations was found in Deo ct al. (2010). 

Hiirvich and Wang (2010, 2009) did not derive limit laws for the log price series or 
the OLS estimator of the cointegrating parameter, but focused instead on properties of 
variances and covariances for log price series and returns, and on lower bounds on the rate 
of convergence for the OLS estimator. 

In this paper, for a slightly modified version of the model of Hurvich and Wang (2010, 
2009), but under assumptions that are more general than theirs, we obtain the limit law 
for the log prices, and for the OLS and tapered estimators of the cointegrating parameter. 
In our result on the limit law for log prices. Theorem 3.1, we allow for a stochastic time- 
varying intensity function in the counting processes. This allows for such effects as dynamic 
intraday seasonality in volatility (as observed, for example, in Deo et al. (2006), as well 
as fixed nontrading intervals such as holidays and overnight periods. We also allow in 
most of our results for asymmetries (leverage effects), and show that this opens up the 
possibility that long memory in durations may affect the rate of convergence of estimators 
of the cointegrating parameter. This raises some heretofore unrecognized ambiguities in 
the choice of a definition of standard cointegration. Finally, we consider the properties 
of the ordinary least squares estimator in a spurious regression, i.e., in the absence of 
cointegration. 

The remainder of this paper is organized as follows. In Section 2 we write the model for 
the log price series and state our assumptions on the counting process, the time-deformation 
functions, and the return shocks. In Section 3, we provide our main results on: the long-run 
behavior of the log-price process (Subsection 3.1), the OLS estimator for the cointegrat- 
ing parameter under weak fractional, strong fractional and standard cointegration (Sub- 
section 3.2), a tapered estimator under weak fractional, strong fractional and standard 
cointegration (Subsection 3.3), a continuously-averaged tapered estimator under strong 
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fractional and standard cointegration (Subsection 3.4) and the ordinary least squares esti- 
mator in the spurious regression case (Subsection 3.5). Section 4 provides proofs. 



2 Transaction level model 

As in Hurvich and Wang (2010, 2009), we consider a bivariate pure-jump transaction- level 
price model that enables cointegration. We define the log-price process y = {yi, ^2) = 
{y{t):t>0)hy 

Niit) A^2{ti, ivi{t)) 

= ^2,A:, (2.1) 

k=l k=l 

N2{t) A^l(t2,iV2{t)) 

y2it)=J2(^^,k + V2,k) + 0-' J2 ei.fc' (2-2) 

k=l k=l 

where for i = 1, 2, Ni{-) are counting processes on the real line (see Daley and Vcrc-,Joncs 
(2003, page 43)) such that, for t > 0, Ni(t) := iVj(0, t] gives the total number of transactions 
of Asset i in (0, t], and ti^k is the clock time (calendar time) for the kth transaction of Asset 
i, with ■ ■ -ti^-i < tifi < < ti^i < ti^2 ■ ■ ■ ■ The quantity N2{ti^Ni{t)) denotes the number of 
transactions of Asset 2 between time and the time ti,Ari(t) of the most recent transaction 
of Asset 1, with an analogous interpretation for Ni{t2,N2{t)) ■ The efficient shock sequences 
{^i,k}'kLi model the permanent component and the microstructure noise sequences {?7j,fc}fcli 
model the transitory component of the log-price process. Efficient shock spillover effects 
are weighted by 6 and 6~^, thus yielding cointegration with cointegrating parameter 6, 
assumed nonzero. A detailed economic justification for this model, derivation of a common- 
components representation, as well as a comparison with certain discrete-time models, is 
given in Hurvich and Wang (2010, 2009). 

Following Daley and Vcrc- .Jones (2003, page 47), a point process is said to be simple if 
the probability is zero that there exists a time t at which more than one event occurs. We 
do not assume that the counting processes are simple. Thus we allow for the possibility 
that several transactions may occur at the same time. The transaction times ti^k are related 
to the point process by the following duality. 

Ni{t) = k^ U^k < t < . 

The durations are then defined for k > 1 hj 

If the process is simple, then Ni{ti^k) = k. Otherwise, it only holds that Ni{ti^k) > k. We 
need the following ergodicity-type assumptions. 
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Assumption 2.1. The sequences {ti,/c} are nondecreasing and there exists Xi G (0, oo) 
such that 

ti^k/k 4 l/A, . (2.3) 

p 

When the counting processes are simple, this is equivalent to Ni{t)/t — ^ A^. Since we 
do not assume simplicity, we must introduce an additional assumption. 

Assumption 2.2. Ni{t)/t — )■ Aj. 

If the counting processes are defined from stationary ergodic durations, tlien Assump- 
tion 2.1 holds. If the couting processes are moreover simple, then Assumption 2.2 also 
holds. Thwo such examples where both Assumptions 2.1 and 2.2 hold are LMSD processes 
and ACD processes, under conditions specified presently. Here and elsewhere, we omit the 
i subscript when this does not cause confusion. 

• For the LMSD process, suppose that = eue"^^ where a is a positive constant, 
{efc} is an iid sequence of almost surely positive random variables with finite mean 
and {Yk} is a stationary standard Gaussian process, independent of {cfc}, such that 
limfc_^oo cov(yo5 ^fc) = 0. It follows from the latter assumption and Guassianity that 
the process [Yk] is ergodic Ibragimov and Rozanov (1978), hence so is {t^}. Thus 
Assumption 2.1 holds with = E[ei]e'^ and since the durations are almost surely 
positive, so does Assumption 2.2. Henceforth, we will further assume that either 
E^i |cov(Fo, Yk)\ <oo or cov(Fo, Yk) ~ ck'^"--'^ for some c> and G (1/2, 1). 

• The ACD model proposed in Engle and Russell (1998) is 

Tk = if^kek, V'fc = w + aTk-i + (3il)k-i, k eZ, (2.4) 

where u > and a, (3 > 0, {^k}'kL-oo is an iid sequence with eo > almost surely, and 
E[eo] = 1. The sequence {Tk}'kL-oQ is then strictly stationary provided such a solution 
to (2.4) exists. From Grimmctt and Stirzaker (2001, Section 9.5), it follows that 
stationary ACD durations are ergodic if E[tq] < oc. The strictly stationary solution 
is then determined hj ipk = Sj^i Ilti (o^^fc-i + and exists if -E[ln(Q;eo + (3)] < 0, 
following arguments given in Bougerol and Picard (1992). As in Aue et al. (2006), 
one can derive that now E[to] is finite if a + /3 < 1 which, on account of Jensen's 
inequality, is also sufficient for E[\n{aeo + /?)]< to hold. Thus, {rk} is ergodic as 
long as a + /3 < 1. In that case. Assumptions 2.1 and 2.2 hold with A = (E[ri])~^ = 
{1-a- l3)/u. 

We now give an example (which is allowed as a special case of our theorems) of how time 
deformations can be used to obtain a nonstationary, possibly non simple point process from 
a stationary simple one. Let Aj(-) be a simple, stationary and ergodic counting process on 
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M with intensity G (0, oo) and let fi be deterministic or random functions such that fi 
is nondecreasing and has cadlag paths with probabihty one. For i = 1,2, define then 

If the functions fi are random, we assume moreover that they are independent of the 
counting processes Ni. 

The functions /j are used to speed up or slow down the trading clock. To incorporate 
dynamic intraday seasonality in volatility, the same time deformation can be used in each 
trading period (of length, say, T), assuming that fi{t) has a periodic derivative (with period 
T and with probability one), for example, fi{t) = t + .5 sin(27rt/T). Fixed nontrading 
intervals, say, t G [Ti,T2), could be accommodated by taking fi(t) = fi(Ti) for t G [Ti,T2) 
so that fi{t) remains constant for t in this interval, and then taking fi(T2) > fi(Ti) so that 
fi(t) jumps upward when trading resumes at time T2. The jump allows for the possibility of 
one or more transactions at time T2, potentially refiecting information from other markets 
or assets that did trade in the period [Ti,T2). 

If the values of series i are only recorded at specific time points (e.g., quarterly in the 
case of certain macroeconomic series) this could be handled by taking fi{t) to be a pure- 
jump function. This provides scope for considering two (or more) series, some of which are 
observed continuously, others at specific times, though not necessarily contemporaneously. 
In future work, we hope to explore this scenario in detail, and its possible connections with 
the MIDAS methodology, see Ghysels et al. (2006). 

The use of the time-varying intensity function fi renders the counting process Ni non- 
stationary. Since it is possible that fi has (upward) jumps, the may also not be simple 
even though the Ni are simple. We now show, however, that Assumptions 2.1 and 2.2 hold 
under the present assumptions. 

Lemma 2.1. Assume that f is a nondecreasing (random) function such that t'^f{t) 
7 G (0, 00) and 

sup|/(t)-/(ni <c 

f>0 

with probability one, where C G (0, 00) is a deterministic constant. Let N be a stationary 
ergodic simple point process, so that Assumptions 2.1 and 2.2 hold for some A > 0. Let N 
be the counting process defined by N{-) = N{f{-)). Then Assumptions 2.1 and 2.2 hold for 
N with A = A7. 

In order to show that our results on estimation of the cointegrating parameter (under 
weak fractional and standard cointegration) hold in this time deformation framework, we 
will in Lemma 4.7 make further assumptions on the fi. These assumptions mathematically 
embody natural economic constraints, viz. minimum duration of trading and nontrading 
periods, maximum duration of nontrading periods and non stoppage of trading time during 
trading periods. 

We now state our assumptions on the return shocks. 
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Assumption 2.3. The efficient shocks {ci^k} o'^e mutually independent i.i.d. sequences 
with zero mean and variance af^. 

Although many of our results would continue to hold if the iid part of Assumption 2.3 
were replaced by a weak- dependence assumption, we maintain the iid assumption here in 
keeping with the economic motivation for the model as provided by Hurvich and Wang 
(2010) that in the absence of the microstructure shocks and in the absence of any depen- 
dence of the efficient shocks on the counting processes, each of the log price series would 
be a martingale with respect to its own past. Since the trades of Asset 1 are not synchro- 
nized in calendar time or in transaction time with those of Asset 2, it seems reasonable to 
assume that the two efficient shock series are mutually independent, as we have done in 
Assumption 2.3. 

The following assumption implies that the microstructure noise does not affect the 
limiting distribution of the log prices. We use =^ to denote weak convergence in the space 
X'([0, oo)) of left-limited, right- continuous (cadlag) functions, endowed with Skorohod's Ji 
topology. See Billingslcy (1968) or Wliitt (2002) for details about weak convergence in 
©([0, oo)). Whenever the limiting process is continuous, this topology can be replaced by 
the topology of uniform convergence on compact sets. 

Assumption 2.4 (Microstructure Noise). The microstructure noise sequences {rji^k} sat- 
isfy n-^''^Y}k=im,k ^ 0. 

Dependence between the counting processes and return shocks allows for leverage effects 
(for example, a correlation between a return in one time period and a squared return in a 
subsequent time period). A transaction- level model yielding a leverage effect was proposed 
(but justified only with simulations) in Hurvich and Wang (2009). Models where the point 
process need not be independent of the return shocks were discussed in Prigent (2001) in 
the context of option pricing with marked point processes. 

We do not make any assumption of independence between the counting processes and 
the microsctucture shocks, unless explicitly noted otherwise. We will, however, assume 
that the counting processes are independent of the efficient shocks except in Theorem 3.1 
and in Section 3.5. 

Assumption 2.4 is all we need to assume about the microstructure noise in order to 
obtain a limit law for the log price series (such as Theorem 3.1 below). However, in order 
to discuss properties of estimators of the cointegrating parameter it is necessary to make 
more specific assumptions on the degree of cointegration. In Hurvich and Wang (2009, 
2010), three different cases were considered, according to the strength of the memory of the 
microstructure noise sequences. These cases were labeled as weak fractional cointegration, 
strong fractional cointegration and standard cointegration. In the current context, where 
there may be a dependence between return shocks and counting processes, special care is 
needed in defining the strong fractional and standard cointegration CclSGS, clS long memory 
in durations may affect the rate of convergence of estimators of the cointegrating parameter 
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in these cases. On the other hand, we will define weak fractional cointegration essentially 
as in Hurvich and Wang (2010). 



Assumption 2.5. The shocks {ei^k}T=-oc> {(^2,k}T=-oc' {Vi,k}T=-oc (^^d {v2,k}T=-oo (^^e 
mutually independent. 

Mutual independence of the efficient and microstructure shock series of a given asset 
can be justified on economic grounds, and is often made in the econometric literature for 
calendar-time models. See, e.g., Barndorff- Nielsen et al. (2008). Mutual independence of 
the two microstructure series is justified by the lack of synchronization of the trading times 
of the two assets. 

We now discuss the weak fractional cointegration case. For H G (0, 1), let Bh denote 
the standard fractional Brownian motion (FBM) with Hurst index H, i.e. the zero mean 
Gaussian process with almost surely continuous sample paths and covariance function 

coy{Bh{s), Buit)) = ]^{s^^ -\t-sr + e"] . 

For H = 1/2, Bi/2 is the standard Brownian motion. 

Assumption 2.6 (Weak Fractional Cointegration). There exists HE (0, 1/2) such that 

k=i 

where Ci,C2 are nonnegative constants, not both zero and B^^^ and B^^^ are independent 
standard fractional Brownian motions with common Hurst index H. 

Under Assumption 2.5, the independence of all the noise series implies that all the 
previous convergences hold jointly. The situation where one of the constants Ci or C2 is 
zero could arise naturally if the memory in one of the microstructure noise series is weaker 
than for the other. 

In the case of weak fractional cointegration, we can define the memory parameter of 
the microstructure noise series as dr, = H — 1/2 G (—1/2,0), and the degree of fractional 
cointegration (i.e. the rate of convergence of partial sums of the cointegrating error) is 
completely determined by c?^. More precisely, in this case the difference between the 
memory parameters of the series of log prices and the cointegrating error (observed, say, 
at equally-spaced intervals of calendar time) yi{j) — Oy2{j), is —drj. This holds regardless 
of any dependence between the counting processes and the microstructure shocks. 

Next we discuss strong fractional and standard cointegration. We start with the as- 
sumption that, for i = 1,2, rji^k = ^i,k — ^i,k-i where {^i,k} satisfy sup^ E[,^j^j!.] < oo. It then 
follows that the cointegrating error at time j is 
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yiU) - Oy2U) 

NiU) N2{j) 

= 5Z ^1'^ ~ ^ 5Z + 6,7ViO) - Cl,0 - 0{^2,N2ij) - ^2,o)- (2.5) 

k=Ni (ta, 0) )+l fc=A^2 (ii,iVi U) )+l 

Under the assumptions we will make in this paper, and also under the assumptions made 
in Hurvich and Wang (2010, 2009), the first two terms on the righthand side of (2.5) are 
weakly dependent, so the degree of cointegration is determined by the rate of convergence 
of partial sums of ^i,Ni{j)- 

Thus we will need to study the sequence ^i,Ni(j)- We do not assume that the microstruc- 
ture shocks are independent of the counting processes. Thus, even if the microstructure 
shocks have zero mean, it may hold that E[^j^Ar-(j)] ^ 0. 

Assumption 2.7 (Strong fractional and standard cointegration). The micro stucture noise 
sequences {rji^k} can he expressed as rji^k = Ci.fc — ii,k-i- There exist H G [1/2, 1) , constants 
^l, ^2 ^'^'^ nonnegative constants Ci,C2, not both zero, such that 

[n-] 
k=l 

where B^^ and B^^ are independent fractional Brownian motions with Hurst index H . 

The case H > 1/2 corresponds to strong fractional cointegration whereas the case 
if = 1/2 corresponds to standard cointegration. 

It might be hard to verify Assumption 2.7 unless the durations are integer valued. Since 
commonly-used duration models do not have integer-valued durations, we will introduce a 
modification of the estimators which involves integrals instead of sums, thus allowing us to 
avoid this restriction. This change requires a corresponding modification of Assumption 2.7. 

Assumption 2.8 (Strong fractional and standard cointegration). The microstucture noise 
sequences {rji^k} can he expressed as rji^k = ii,k — ^i,k-i- There exist H G [1/2, 1), constants 
/ig and nonnegative constants Ci,C2, not hoth zero, such that 

pn- 

Jo 

In their strong fractional cointegration case, Hurvich and Wang (2010) assumed, for 
G (—1,-1/2), that cov {^i^k,^i,k+j) ~ Kj^'^^~^ as j' — )■ oo where K > and = 
+ 1 G (0,1/2). They then showed (in their Lemma 3), under the assumptions made 
there, that cov(,^j jv.(fc), 7Vi(fc+j)) ~ as j — oo where K' > 0, so that the degree of 

fractional cointegration was completely determined by the rate of decay of coy {^i^k, ^i,k+j)- 
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However, the proof of this result rehed on the assumption that the microstructure shocks 
are independent of the counting processes, an assumption which we do not make here. 

We next provide an example showing that under dependence between the microstruc- 
ture shocks and the counting processes, it is possible for {^i^k} to be weakly dependent, and 
yet the rate of convergence of suitably normalzed integrals of the process {^i,Ni(t) '■ t > 0) 
is determined by the degree of long memory in durations. Suppressing the i subscript, we 
have the following lemma. 

Lemma 2.2. Suppose that {rk} is given by the LMSD model Tk = ^k^'^ , {cfc} i.i.d. 
standard exponential, independent of the stationary standard Gaussian series {Yk}, which 
satisfies cov(yo,^n) ~ cn"^^^'^ where c > 0, and G (1/2,3/4). Define = Y^_^_^ — 1. 
Then the autocovariance function of {^k} is summable and there exists c' > such that 

n-''^Y.^k^cB . (2.6) 

k=l 

Nevertheless, the randomly-indexed continuous-time process ^N{t) has long memory in the 
sense that there exists a constant fi* , as well as a constant c" , such that 

n-""^ / Hn^s) - ds c"Bh^ . (2.7) 

Lemma 2.2 shows that long memory in durations can induce the same degree of long 
memory in the cointegrating error (2.5) in calendar-time, even though the microstructure 
shocks, which are the source of the cointegration, have short memory as a sequence in 
transaction time. In Lemma 2.2, this phenomenon was achieved by imposing a particular 
functional relationship between the (zero mean) microstructure noise and the persistent 
component of the durations, = ^k+i ~ This relationship implies a leverage effect, 
since coii[^kiTk+i) = 1 / A/2(e — 1) ~ .539 > 0. In other words, a strongly negative mi- 
crostructure shock to the return leads to a tendency of the next observed duration, as well 
as subsequent durations, to be shorter than average. Such a string of short durations in- 
creases the volatility, e.g., the expectation of squared calendar-time returns, as shown, for 
example, under a particular return model in Deo et al. (2009). Furthermore, evidence that 
stock market intertrade durations have long memory was provided in Deo et al. (2010). 

In the absence of dependence between the counting processes and microstructure noise 
series, in both cases of strong fractional and standard cointegration, the memory of du- 
rations cannot affect the memory of the cointegrating error. See Lemma 4.9 for strong 
fractional cointegration and Lemma 4.10 for standard cointegration. 
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3 Main results 



3.1 The long-run behavior of the bivariate log-price process 

With the assumptions made in Section 2, the long-run behavior of the bivariate process 
y = {yiiy2) can be determined. The following theorem shows that the log-prices are 
approximately integrated. Even though independence is assumed between the various 
shock series, the log-price process y = {y(t): t > 0) exhibits a nontrivial variance-covariance 
structure which is determined by a complex interplay of the model parameters. 

Theorem 3.1. Under Assumptions 2.1, 2.2, 2.3, 2.4, n~^^'^{yi{n-),y2{n-)) =^M, where 

B = (ai,eV^fil + 9(r2,eV^2B2 , ^"Vi,eV^fil + a2,eV^B2) . (3.1) 

and Bi and B2 are independent standard Brownian motions. 

In Theorem 3.1, we have not assumed that the counting processes are independent of 
either the efficient shocks or the microstructure shocks. 

Hurvich and Wang (2010, 2009) have in their Theorem 1 computed the long-run vari- 
ances of yi(t) and ?/2(t) which are given as {(jf^Xi + 9^(T2f,\2)t and {af^\i/9^ + cr|gA2)t, 
respectively. Our theorem yields the variances as well as the covariances in the limiting 
distribution of {t~^^'^y{t) : t > 0). More importantly, our theorem provides the limiting 
distribution itself for the (normalized) log-price process y which, in turn, can be used for 
asymptotic statistical inference. Indeed, most of the subsequent results in this paper use 
Theorem 3.1 and its proof as a building block. In particular, a slightly generalized ver- 
sion of this theorem directly yields asymptotics for estimators in spurious regressions and 
therefore can be used to motivate tests for the null hypothesis of no cointegration, as we 
describe in Section 3.5. 

3.2 OLS estimator of the cointegrating parameter 

In this section, we derive the asymptotic behavior of the ordinary least-squares estimator 
(OLS) of the cointegrating parameter 6. To do so, we assume that the log-price series are 
observed at integer multiples of At. We will work here, without loss of generality, with 
At = 1 in order to keep the notation simple. Then (2.1) and (2.2) become 

Nl{j) '^2{h.N^U)) 

fc=i fc=i 

N2ij) A'lfe.iVaO)) 

y2{j)= ^ie2,k + V2,k) + 0~^ ^ ei,k- 
k=l k=l 
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Regressing yi{l), ■ ■ ■ ,yi{n) on 1/2(1), • • • ,y2{n) without intercept, we obtain the OLS esti- 
mator of 6 as 

Hurvich and Wang (2010, 2009) have shown in their Theorem 6 (under conditions that are 
for the most part stronger than the ones we assume here) that is weakly consistent 
for 6 and obtained a lower bound on the rate of convergence in the case of weak fractional, 
strong fractional and standard cointegration. The exact limit distributions, however, were 
not given. We fill in this gap next for weak fractional cointegration. 

3.2.1 OLS under weak fractional cointegration 

Theorem 3.2. Let Assumptions 2.1, 2.2, 2.3, 2.5 and 2.6 hold. Assume in addition that 
the counting processes Ni and N2 are mutually independent and independent of the efficient 
shocks and there exists a constant C such that 

supE[ti_7v,(s)+i - s\<C . (3.3) 

Then 

jlB{t)BH{t)d.t 



jlB\t)dt 



where B is a standard Brownian motion, Bh is a fractional Brownian motion, independent 
of B , and 

„2\2H I ^2\2H 
y2 ^ ^1^1 + "^2^2 

' e^'x^ai + A24, • 

The result in Theorem 3.2 is similar to that obtained in Robinson and Marinucci (2001, 
Proposition 6.5, formula (6.8)), under their Assumption 6.1, for which a sufficient condition 
(their formula (6.5)) was verified in Marinucci and Robinson (2000) to hold for weak (but 
not strong) fractional cointegration in the case where the process is linear with respect to 
iid innovations. 

Remark 3.1. For the LMSD model described in the previous section, it is proved in 
Lemma 4.8 that (3.3) holds under the assumptions of Lemma 2.2 that pertain to the 
durations. For the ACD model described in the previous section. Lemma 4.6 shows 
that (3.3) holds as long as IE[r|] < 00. Carrasco and Chen (2002, Corollary 6) prove 
that if E[{/3 + aetf] < 1 then E[r|] < 00. For the LMSD or^ACD models, or for any 
duration model such that the corresponding counting process satisfies (3.3), the time 
deformed process A^ also satisfies (3.3), provided the time deformation function / satisfes 
the assumptions of Lemma 4.7. Therefore Theorem 3.2 holds for LMSD or ACD durations 
with deformation functions satisfying the assumptions of Lemmas 2.1 and 4.7. 
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3.2.2 OLS under strong fractional and standard cointegration 



We now consider tlie case wliere tlie microstructure noise series {rji^k} are differences of 
strongly or weakly dependent processes 

Theorem 3.3. Let Assumptions 2.1, 2.2, 2.3, 2.5, and 2.7 hold. Assume moreover that 

• the efficient shocks are i.i.d. Gaussian, 

• the counting processes Ni and N2 are independent of each other and independent 
of the microstructure noise sequences and of the efficient shocks and there exists a 
constant C such that 

supE[(t,,^^(t)+i-t)2] <C. (3.4) 
t>o 

• E[{,,fc] = 0, supfcE[{2J < ^^^^ _ 
Then, 

• ifl/2<H< 1, 



where B is standard Brownian motion independent of the standard fractional Brow- 
nian motion Bh; 

• ifH= 1/2, n(^°^^ -6) = Op{l) . 

The rate of convergence obtained in the standard cointegration case improves on the 
one obtained by Hurvich and Wang (2010). 

The assumptions in Theorem 3.3 are quite strong, ruling out leverage effects and pro- 
viding one motivation for our subsequent consideration of tapered estimators. 

Remark 3.2. For the LMSD model described in the previous section, it is proved in 
Lemma 4.8 that (3.4) holds under the assumptions of Lemma 2.2 pertaining to the du- 
rations. Similarly, for the ACD durations. Lemma 4.6 and Carrasco and Chen (2002, 
Corollary 6) imply that if E[(/3 + aei)^] < 00 then (3.4) holds. However, we are unable to 
show that Assumption 2.7 holds except in certain cases of integer- valued durations. We 
do not pursue this further here. 
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3.3 A Tapered Estimator of the Cointegrating parameter 

Even in existing discrete-time models for cointegration the OLS estimator lacks any par- 
ticular optimality properties. Here we consider an estimator based on discrete Fourier 
transforms of the tapered differences of yi{j), y2{j), 1 < J < It was shown in 
Chen and Hurvich (2003a) that this estimator can have a faster rate of convergence than 
OLS in certain cases of fractional cointegration. In the weak fractional cointegration case, 
our limit results for the tapered estimator (Theorem 3.4) are obtained under identical 
conditions as those assumed in Theorem 3.2 for OLS. However, under strong fractional 
and standard cointegration, the conditions for our results on the tapered estimator (The- 
orem 3.5) allow for leverage, unlike the corresponding theorem for OLS. 

We introduce all relevant notation using a generic time series {xj}'^_^. Let /z: / — t- R 
be a general continuous taper function on an open interval / containing [0, 1] such that 
h{0) = h{l) = 0. For £ = 1,2, ... , denote by ue = 27ri/n the Fourier frequencies. The 
tapered DFT of {xj}^?^^^ with taper function h is defined by 

n . n . 

j=l 3=1 

where he{t) = /i(t)e^'^^^*. Denote by {Axj}'^_^ the first difference of the series {xj}, where 
Axj = Xj — Xj^i. We define the tapered DFT of the first difference {Axj}'jL_^ with taper 
function h by 

n . n . 

dA., = J2h{i)^^^^''^' = Y.^^{i)^'^r (3.6) 

j=i i=i 

In our setting, we observe the cointegrated component processes yi and y2 at equidis- 
tant sample points. Defining the cointegrating error zj = yij — 6y2j and following 
Chen and Hurvich (2003b), we can now introduce the estimator 

e^^P = Re(4) , 

where Iie{z) signifies the real part of a complex number z = a + ib and, letting z = a — ib 
be the complex conjugate of z, 

a _ Y^7=i dAyi,e dAy2,e 

Y^"^ M 12 

2^e=i \ "'Ay2,e\ 

Therein, any tapered DFT of differenced sequences is defined according to (3.6). Note that 
6'^°'^ is the real part of the ratio of the averaged tapered cross-periodogram between the 
series yi and y2 and the averaged tapered periodogram of the series y2- 
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3.3.1 Discrete tapered estimator under weak fractional cointegration 



Theorem 3.4. Let Assumptions 2.1, 2.2, 2.3, 2.5 and 2.6 hold. Assume moreover that 
the counting processes are mutually independent and independent of the efficient shocks 
and (3.3) holds. Then 



1/2-H ^QTap _f^^d^ / cf Af^ + e^clXf' Er=i Re [!l Ms) ^B{s) jl hejt) dBnit) 



e-^Kal,\, + X^al^ ^™ ^ ^^1 ^^^^^ ^^^^ 



where B is a standard Brownian motion, Bh is a standard fractional Brownian motion 
and B and Bh are independent. 

Since the assumptions of Theorem 3.4 are the same as in Theorem 3.2, Remark 3.1 also 
applies here. 

3.3.2 Discrete tapered estimator under strong fractional and standard coin- 
tegration 

Theorem 3.5. Let Assumptions 2.1, 2.2, 2.3 2.5 and 2.1 hold. Assume morover that 
the counting processes are mutually independent and independent of the efficient shocks 
and (3.4) holds. 

• If 1/2 < H <1, then 

^2 I Q2^2 Re ( hAs) dBis) h'As) dB^is^ 



where B^ is a standard fractional Brownian motion independent of the standard 
Brownian motion B. 

• If H= 1/2, n{9l'^^-e) = 0p{l). 

The assumptions of this theorem are weaker than those of Theorem 3.3 on the OLS 
estimator. The microstructure shocks are not assumed to be independent of the count- 
ing processes and the efficient shocks are not assumed to be Gaussian. Theorem 3.3 can 
presumably be proved without the Gaussian assumption. It might be much more diffi- 
cult in the proof of Theorem 3.3 to avoid the assumption of independence between the 
microstructure shocks and the counting processes. 

Remark 3.2 remains relevant here. In particular, we are currently unable to show that 
Assumption 2.7 holds except in certain cases of integer- valued durations. This motivates 
our consideration of a continuous-time tapered estimator in the following section. 
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3.4 A Continuous-Time Tapered Estimator 

The estimators of 9 we have considered so far are based on equally-spaced observations 
of the log price series. However, under the model (2.1), (2.2), a continuous-time record 
is available, and it is of interest to consider using all of the available data to estimate 9. 
Here, for the sake of theoretical tractability, we consider a tapered estimator 9 based on 
continuously-averaged log prices on adjacent non-overlapping time intervals. Since the 
problems with discretization appear only in the strong fractional and standard cointegra- 
tion cases, we only consider them in this section. There is no difference in the case of weak 
fractional cointegration. 

We first establish some notation. Let {X{t)} be any time series defined for all t > 0, 
and suppose that we have data on {Xt} for t G [0,T]. Let 5 > be fixed. In practice, we 
might take 5 to be 5 minutes, but the choice of 5 does not affect the asymptotic distribution 
we derive below. Define n = lT/6\, X{0) = 0, and 

X{k) = / X{u)du, A; = I,-- - ,n. 

Ju={k-1)5 

Then we can define an estimator 9s based on these averaged observations by 



with 



9, 



n,5 



J2T=1 \dAy2,- 



|2 



3.4.1 Continuous-time tapered estimator under strong fractional and stan- 
dard cointegration 

Theorem 3.6. Let Assumptions 2.1, 2.2, 2.3, 2.5 and 2.8 hold. Assume morover that 
the counting processes are mutually independent and independent of the efficient shocks 
and (3.4) holds. 

• Ifl/2<H < 1, then 

_^ , S^^icl + 9^cl) TJLi Re (/o^ ht{s) dB{s) h'^{s)dBH{s] 



n^/-^(Cj ^ ^^^^ I ^, ^^^^^ ^^^^^ 

where Bh is a standard FBM independent of the standard Brownian motion B. 
If H = 1/2, n{9lf-9)=Op{l). 
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Because Assumption 2.8 involves an integral rather than a sum, we are able to verify 
that it holds for certain models with noninteger durations such as ACD and LMSD under 
certain relationships with the microstructure shocks. 

In Theorem 3.6, we allow for leverage effects, and therefore care is required in defining 
standard cointegration. As demonstrated in Lemma 2.2 (which assumes LMSD durations) 
if there is a leverage effect, even when the microstructure shocks are the differences of a 
weakly- dependent sequence, the cointegrating error need not be /(O). In such a case we 
have strong fractional cointegration rather than the standard cointegration which might 
have been anticipated. 

It is also possible that even though a leverage effect exists, the memory of durations 
has no effect on the degree of cointegration. Specifically, if in Lemma 2.2 we replace 
= Y,\, - 1 by ^fc = i/2(n+i) - .75if3(n+i), where H,{y) = - 1 and H,{y) =y^-3y 
(the second and third Hermite polynomials, respectively), then there is a leverage effect 
with corr(rfc+i, ,^fc) = .082. Nevertheless it follows from an argument similar to the proof 
of Lemma 2.2 that Assumption 2.8 holds in this example with H = 1/2, so that we have 
standard cointegration and Theorem 3.6 holds with H = 1/2. 

Lemma 4.10 provides an example of standard cointegration allowing for both time 
deformation and dependence between the counting processes and microstructure shocks. 
Theorem 3.6 would hold for this example with H = 1/2. 

3.5 Spurious Regressions 

In this subsection only, we consider a non-cointegrated version of the model defined by 
(2.1) and (2.2), 

Ni(t) A^2(ii,jvi(t)) 

2/1 (^) = X^(ei,fc + ?7i,fc) +6*21 ^ e2,fc, (3.7) 

fc=l k=l 
N2{t) ^l{t2,N2(t)) 

y2it) = ^{e2,k + V2,k) + di2 ^ ei,fc, (3.8) 

k=l k=l 

where 9i2 ^ ■ We examine here the properties of the OLS estimator in the (spurious) 
regression of yi on y2 in discrete time and then briefly discuss corresponding tests for the 
null hypothesis of cointegration. Corollary 3.1 below follows directly from the proof of 
Theorem 3.1. 

Corollary 3.1. If Assumptions 2.1, 2.2, 2.3 and 2.4 are satisfied and y = {yi,y2) is given 
by (3.7) and (3.8) with 9i2 7^ 621, then as n ^ 00, 

^y{nu) : n G [0, 1]^ ^ By = {By{u) : n G [0, 1]), 
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where By is a bivariate Brownian motion with 2x2 covariance matrix S = (Sjj : i,j = 1, 2) 
given by the entries 



2 



^1,1 - 

Sl,2 = 6'l2Ai CT^ g + 6'2lA2 e = S2,l. 

Next, we consider the discretization of yi(t) and y2(t) given by (3.7) and (3.8) at integer 
time values, 



Ni{j) 7V2(ii,iVi(,)) 

yi,j = ^ (ei,fc + Vi,k) + O21 ^ e2,fc, (3.9) 

k=l k=l 
N2{j) Afl(t2,iV2(j)) 

?/2,j = ^ (e2,fc + ?72,fc) + 6^12 ^ ei,fc. (3.10) 

k=l k=l 

Regressing yi^i, . . . , yi^^ on 2/2,1, • • • , 2/2,n without intercept, we obtain the OLS estimator 



En 
j=i y2,jy i,j 



Corollary 3.2 below follows directly from Corollary 3.1 and the Continuous Mapping The- 
orem. 

Corollary 3.2. If Assumptions 2.1, 2.2, 2.3 and 2.4 are satisfied and y = (?/i,?/2) is given 
by (3.7) and (3.8) with 612 7^ 621, then as n ^ 00, 

. ^ B2,yiu)Bi^y{u) du 
lo Blyiu) du 

where By = (Bi^y, B2^y) is the bivariate Brownian motion given in Corollary 3.1. 



Corollary 3.2 together with Corollary 3.1 can be used to motivate tests for the null 
hypothesis of no cointegration. We do not pursue the details here, but it seems clear 
that the null distribution for unit root tests based on the residuals {yij — 6ny2,j}^=i can 
be derived from Corollaries 3.1 and 3.2, and that these null distributions will have form 
similar to the distributions listed, for example, in Hamilton (1994, Proposition 19.4). 



4 Proofs 

Proof of Lemma 2.1. Since N is stationary and ergodic, there exists A G (,00) such that 
N{t)/t 4 A a St ^ 00, a. s. See Daley and Vere- Jones (2003, diplays (12.2.3) and (12.2.4)). 
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Now, 

Njt) _ Njm) _ Nifit)) m 
t t fit) t 

Since by assumption f{t) oo (in probability if / is random), the assumption implies 
that the first term on the righthand side converges to A. By assumption, the second term 
satisfies t~^f{t) 7. Thus, A, so that Assumption 2.2 holds with A = A7. We 

note that N(t^) < k, thus 

^ ^ Njh) _ ^ ^ Njt,) -k ^ Njtk) - Njq) 
~ k k ~ k 

_ ^ ^ NjfiU)) - iV(/fa-)) ^ ^ ^ Njfih)) - NjfiU) - C) 
k ~ k 

using the definition of N and the boundedness requirement on the jumps of /. Since the 
process is stationary, taking expectations, we have 

E[iV(/(t,)) - Nifih) - C)] = E[iV((0, C])] < 00 . 

(In the case where / is random and independent of this holds by taking conditional 
expectation first). Thus it follows that 

^ = 1+0.(1/.). 

Now, this implies that N{tk)/tk converges in probability to 1. Thus 

h ^ tk Njtk) P 1 
k N{tk) k ^ A ■ 

□ 



Proof of Lemma 2.2. Denote Hr by H to simplify the notation. Define C,k = ^k+i ~ 
1. Then C,k is centered, has finite variance summable autocovariance function, since 
cov{^o,^k) = 2cov^(yo, ^fc+i)- Thus {^k} has a summable autocovariance function because 
H e (1/2,3/4). By Arconcs (1994, Theorem 4), this implies that {^k} is in the domain of 
attraction of the standard Brownian motion, i.e. 

[n.] 
k=l 

with c'^ = var(^o) + 2 ZlfcLi cov(^o, ^k)- This proves (2.6). 

Assume now that = e/je'^^*^ (with a = 1 in the statement of the Lemma). The 
properties of Hermite polynomials yield that ¥.[e"^° Hj(Yo)] = a^e'^ for all j > 1. Denote 



19 







now = E[rk] = E[e'^^^] = e'^'/^ m = E[^k-irk] = E[(Fj2 _ 1)6^.] = ^rV'/^ and 
Cr{y) = {y"^ — l)e'^^ — m. We now prove that (2.7) holds with /i* = Am. Write 

T 

(^7V(.) - Am) ds 

Af(T) 

= ^ Tk+l^k - y^mT + (tAr(r)+l - T)^N{T)+1 
N{T) N{T) 

= J2 (^fc+i - ^)^ke''^''^' + G{Yk+i) + m(iV(r) - AT) - (tjv(T)+i - T)^NiT)-,i • 

fc=0 fc=0 

(4.1) 

By Lemma 4.8 and applying Holder's inequality, it can be shown that (t iy(^T)+i—T) {C,n(t)+i — 
p) = Op(l). Since the sequence {efc} is independent of the Gaussian process {Y^}, the first 
term in the righthand side of (4.1) is in the domain of attraction of the standard Brownian 
motion, and the normalizing sequence is ^/n. Thus we must obtain the joint asymptotic 
behaviour of Efif^ G{Yk) and N{Tt) - XTt. 

The durations are in the domain of attraction of the fractional Brownian motion with 
Hurst index H, since 

E(r. - A-^) = Y^i^, - l)e-^^ + X:(e-^^ - A^) . 

k=l k=l k=l 

The first term in the righthand side is Op{^/n) and the second sum, suitably normalized 
converges to the fractional Brownian motion with Hurst index H because the function 
X — >■ e'^^ — has Hermite rank 1. See e.g. Arcones (1994). More precisely, let ci = 
E[Yie'^^^] = cre'^ and define g{y) = e'^^ — A~^ — ciy. The function g has Hermite rank 2, 
and since H G (1/2,3/4), this implies that 



var 5^^7(n) =0{n) 



Thus X^fc=i(Tfc — A~^) is asymptotically equivalent to Ci Yl^^i Y^. Let Bu denote the stan- 
dard fractional Brownian motion with hurst index H. The assumption on the covariance 
of the Gaussian process {yu} implies that 

[n.] 
k=l 

with if^ = c/{H{2H - 1)}. Denote now ca = E[YiG{Yi)] = a{a^ + 2)e'^'/2 ^^^^ ^gg^g 
h{y) = G{y) — C2y- Then h has Hermite rank 2 and thus by similar arguments as above. 
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Ylk=i ^(Yk) is asymptotically equivalent to C2 Ylk=i Thus we obtain 



n 



[nt] [nt] 



k=l k=l 



By Vervaat's Lemma (see Vervaat (1972) or Resnick (2007, Proposition 3.3)), the previous 
convergence implies that 



n 



[nt] 

I Mint)- Xnt, J2G{Yk) 1 {-\ci^BH{\t),C2^BH{t)) . 



k=l 

By the continuity of the composition map, this yields 



N{nt) 

N{nt) - Xnt, GiYk) 1 ^ (-XciifBHiXt), c^^BniXt)) . 



k=l 

Next we obtain that 



n 



{N{nt) 
J2 G{Yk) + m{N{nt) - Xnt) ^ ^{c2 - Xmci)BH{Xt) 
k=l 



with C2 — Xmci = 2cre'^^/^ > 0. We conclude that n ^ {^Nis) — Xm} ds =^ (f{c2 — 
Xmci)BH- □ 



4.1 Proof of Theorem 3.1 and Corollary 3.1 

We first need the following Lemma. 

Lemma 4.1. Under Assumption 2.1 and 2.2, Ni{tj^j:q.(nt))/n converges in probability uni- 
formly on compact sets to Xit, where = {1,2}. 

Proof of Lemma 4-1- The sequence of (random) functions Ni{n-)/n is nondecreasing and 
converges pointwise in probability to Xit by ergodicity. A sequence of nondecreasing func- 
tion converging to a continuous function converges uniformly on compact sets. This results 

is known as Dini's Theorem. Cf. Resnick (1987, page 3). Thus the convergence of Ni{n-)/n 

p p 

is uniform on compact sets. Assumptions 2.1 and 2.2 imply that Ni{t) — t- oo and tj^„ — t- oo. 
Thus 

Njitj^Njinu)) _ Niitj^N^jnu)) tj^N.jnu) Njjnu) P 1 ^ . _ . 

n h,Nj{nu) Nj{nu) n Xj 

Applying again Dini's lemma, we also have that Ni[tj^fyr.(^nu))/n converges uniformly on 
compact sets to XiU. □ 
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Proof of Theorem 3.1. Denote S^ .^{t) = Y^^k=i^hk and S^,^{t) = Y^^k=iVi,k, i = 1,2. Under 
Assumptions 2.3 and 2.4, n~^^'^{Sln, 5*1 S^^, S^j^) converges weakly to {ai^eBi, cr2,e-B2, 0, 0), 
where Bi and B2 are independent standard Brownian motions. This follows from the in- 
dependence of ei and 62 and the local uniform convergence to in probability of n~^^'^Sl'^. 
With the previous notation, (3.7) and (3.8) become 

yiint) = Sl,iN,{nt)) + e2i^2%(A^2(ti,iv,(n*))) + Sl^{N,{nt)) , 
y,{nt) = Sl^{N2{nt)) + euSl^{N,{h,N,int))) + Sl^{N,{nt)) . 

By Lemma 4.1 and the continuity of the composition map on C x C endowed with the 
metric of uniform convergence on compact sets (see e.g. Billingsley (1968, Chapter 3, 
Section 17)), we obtain the joint convergence of 

n-'/' {Sl^{N,in-)),SUN,{h,N,in.))), 

SUmn-)), 52%(iV2(ti,^,(„.))), SlXN,{n.)), Sl^iN.in.))) 

to (ai^eV%"-Bi5 f^i,eV^-Bi5 f^2,eV^-B2, cr2,eV^-B2, 0, 0). This yields Corollary 3.1 and The- 
orem 3.1 by setting 6*21 = 6 and 612 = 0~^. □ 



4.2 Proof of Theorems 3.2 and 3.3 

Proof of Theorem 3.2. Write 

2 E?=i{yi(j) -^y2(j)}y2(j) 

Ei=i?/2U) 

Assumptions 2.1, 2.2, 2.3, 2.5 and 2.6 imply those of Theorem 3.1. Thus we can apply the 
Continuous Mapping Theorem and obtain 

n „1 

r^-'E2/2(j)^{^"%<e + A24e} / B\s) As , (4.2) 
3 = 1 

where i? is a standard Brownian motion. Thus, in order to study the convergence of ^„ — ^ 
suitably renormalized, it suffices to study the sum 

n 

E{yi(j)-^y2(j)}y2(j)- 
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We further decompose the cointegrating error. Denote 

Niij) ^^2(^l,JVlO■)) Niij) 

fc=l k=l k=l 

NiU) A^i(i2,iV2 0)) N2{j) 

VlU) = Y ^2,fc + ^"^ Y ^1'*^ ' ^2(i) = Y ^2,fc ' 
fc=l fc=l fe=l 

Ni(j) N2U) 
fc=A^i{t2,iV2 0-))+l fc=A^2{tl,iViO))+l 

With this notation, we can write 

n n n 

- ^2/20')}y2(j) = Y^n, - er,,}y,{j) + Yivlij) - %2^(j)}?/2(j) . (4.3) 
j=i j=i j=i 

Applying Theorem 3.1, Assumption 2.6 and the Continuous Mapping Theorem, we obtain 

n 

^"'~''Y.^ym-dyl{j)}y2{j) 



n 



A [ {rlv^ai,e5i(t) + V^a2,eB2(t)}{Ci5i,H(Ait)-^C2B2,H(A2t)}dt 

= S / B{t)BH{t)&t 
Jo 

where 5 is a standard Brownian motion, Bh is a fractional Brownian motion, independent 
of B and 

= (r^Aia,% + A24e)(c?A?^ + ^^c^Af ) . (4.4) 
There only remains to prove that, for i = 1,2, 

n 

"-'/'E^M^2(j) = 0p(l). (4.5) 

i=i 

The convergence of n"^^'^y2 is uniform on [0, 1], so n"^/^ maxi<j<„ |y2(j)| = Op{l). There- 
fore, it suffices to prove that 

n 

^''EK.I=0p(1). (4.6) 
Recall that Ni{s) < k 4^ U^k > s. Thus, for k < Ni{n), 

Ni{t2,N2(j)) <k< Ni{j) ^ t2,N2{j) < tl,k < j ■ 
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The first inequality on the righthand side means that there is no point of between ti^k 
and j, i.e. j < t2^N2(ti Let A2{t) = t2^N2it)+i ~ t denote the forward recurrence time of 
N2, i.e. the time between t and the next event of after t. Thus, 

n n N-i(j) Ni{n) 

^Yl M = Yl M{Mti,k) + 1} . 

j=l j=l fc=Afi{t2,jv2 (,))+! k=l 

We thus get the bound for the conditional expectation given the sigma-field J\f generated 
by the counting processes A''i and N2: 

Ni{n) 



E 



Conditioning on A''i and applying (3.3) yields 



k=l 



E 



< CNi{n) = Op{n) 



This proves (4.6) and concludes the proof of Theorem 3.2. 



□ 



Proof of Theorem 3.3. The proof is a consequence of the convergence (4.2), the decompo- 
sition (4.3), and Lemmas 4.2 and 4.3, whose assumptions are those of the Theorem. □ 

Lemma 4.2. Under the assumptions of Theorem 3.3, 



n 



n „i 
3 = 1 



(4.7) 



where Bh is a standard fractional Brownian motion independent of B and 

j:o = ie-'Xial + X2al,){cl + e'cl). 

Proof of Lemma 4.2. Denote Sn = Ej=i{2/i (j) " ^^2 (j)}2/2(j) and write y2 = yl + with 
obvious notation. Denote Q = Uiij) — Oy2{j) = C,i,Ni{j) — ^^2,N2ij)- Then 

n n 

5n = 5^0l/2(j) + E0W). (4.8) 

i=i i=i 

By the last part of Assumption 2.7, the last term in the righthand side of (4.8) is Op{n). 
Consider the first term in the righthand side of (4.8), say Si^n- Write 

n N2U) n ^i{*2,iV2U)) 

j=l k=l 3=1 k=l 

N2{n) ^l{t2,N2(n)) 

= E E + E ^1.^ E 

k=l {j<n:N2ij)>k} k=l {i<n:Afi(t2,]V2(j))>fc} 

= Tin + ^T2n . 
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Denote Wn{t) = Yl^^i Cj- Since N2{j) < k iS j < t2,fe, we obtain 



N2{n) 



k=l 



By Assumption 2.7 and Theorem 3.1, n^^/^^^y^''{n)Wn{l) A v^(72fi2(l)^(l) with Z = 
CiB^^^ — 6c2B^^^ = cl + Bh- Let the last term be denoted by Un- Since the shocks 
Ci^k are i.i.d. Gaussian, we can compute the characteristic function of Un- 

E[exp{itn"^/2~^f/„}] = E 



E 



exp 



exp 



-.2 j-2 -I ^2{n) 

1 ^ (^-^^^(,^^^/,))^ 



2 n 

fc=i 



Z^(s)ds 







The convergence is actually joint with that of n ^y2^{n)Wn, thus we have 

n-^/'-^Ti,„ ^ V^(72,efi2(l)^(l) - V^^2,e / ^(S)d52(s) . 

Jo 

The limit can also be written as V^c'"2,e /q^ B2{s) dZ{s). Consider now the term T2^„. Note 
that Ni{t2,N2{j)) <kiSj< t2,7V2(ti,fc)+i- Thus 

^l(*2,JV2(n)) ^l(*2,JV2(n)) 

T2,n = ^ ei,fciyn(l) - ^ ei,fcW'„(t2,Ar2(ti,fe)+l/^) • 
fc=l A:=l 

By similar arguments as previously, we obtain 



n 



''"^'^T2,n ^ V^^lfil(l)Z(l) - a/aT^I Z{S) d5i(s 

^0 



All convergences hold jointly, thus (4.7) holds. 
Lemma 4.3. Under the assumptions of Theorem 3.3, 

n 

^{rij - Qr2,j)y2{3) = Op{n) , 
Proof of Lemma 4-3. We first study the term with rij and split it into three parts. 

n n n n 

Ynjy2{j) = ^njy^'ij) + Yrijy^^{j) + ^rij2/^(j) 
j=i j=i j=i j=i 



□ 



(4.9) 
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We start with the last one. Recall that A''i(t2,Af2(i)) < k < Ni{j) iff ti^k < j < ^i,fc + ^2(^i,fc)- 
Thus 



i=i i=i A^i(t2,]V20))<fc<^i(i) 

If the microstructure shocks are independent of the counting processes, then 



(4.10) 



E 



A'l(n) 

fe=l 



X ^2,7V2(i) I I A/" 

V<l,fe<i<<l,fe+A2(tl,fc) 



A'i(n) 



<CX(^2(ti,fc) + l)'supEg,]. 
Conditioning on A^^i and then taking expectation yields 

^[(E^i-'^ X ^2,JV20)) ] <CE[iVi(n)]supE[{l + A2(t)}2]supE[e,2] = 0(n). 

J=l ti.k<j<ti.k+A2{ti,k) * ^ 

Consider now i?2,n = Z]"=i ^i,il/2'(j)- 

j=l A'i(i2,iV2U))+l '==1 il,fe<J<il,fe+A2(ii,fe) 

By independence of the efficient shocks and the counting processes, we have 



Ni{n) 



E[4„ I AT] < CN,{n) J2 iMti,k) + 1)' = Op( 



n 



k=l 



This proves that -R2,n = Op{n). Consider finally Rin = J2]=i^i,jy2^U)- ^y definition, 
ei^fc is independent of 1/2^ (j) for j such that A''i(t2,7V2 0)) < k. Thus, we can compute the 
conditional variance given A/". 



Ni{n) 

E[Rl^ I Af] = X ^ 

k=l 



Jl,k<j<tl,k+A2{tl,k) 



Nlin) 



< CN2{n) E {M{Kk) + 1)' = Op[ 



n] 



k=l 



by (3.4). This concludes the proof of Lemma 4.3. 



□ 
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4.3 Proof of Theorems 3.4 and 3.5 



Write 



(^n — t7 H ^^Tfi — r-; tt; h 



with = yl{ji) - 6'2/^(j), r(j) = ri(j) - 6'r2(j) and 
By summation by parts, since h(^) = h{l) = 0, for any time series {xj}, we can write 



Ni(j) N2{j) 

ei,k , r2{i) = ^ e2,k ■ 

fc=A'lfe,]V2(j)) + l fc=A'2(tl,]Vi(j)) + l 



n-1 



n-1 



c^Ax/ = X]{/iKj7'^) - hi{{j + l)/n)}xj = ^Wi{j,n)xj (4.11) 

j=0 ^ 



j=0 



with we{j,n) = n{hg{{j + l)/n) — he{j/n)}. Applying (4.11) to ?/2 yields 



dAy2,e = Y] we{j, n)y2{j) 



Since the assumptions of Theorems 3.4 and 3.5 imply those of Theorem 3.1, the Continuous 
Mapping Theorem yields 

{n-^^^dAy2/, 1 < ^ < A ^ Kis)B{s) ds , 1 < ^ < m| (4.12) 

where i? is a standard Brownian motion and = 9~'^\ia\^ + \2<J2e- -^Y integration by 
parts, the integral can also be expressed as 

( h'^s)B{s)ds= [ hi{s)dB{s) . 
Jo 



This in turn implies 



n 



£=1 e=i -^^ 



(4.13) 



Applying now (4.11) to we obtain 

^ n—l 

dAyv,e = — y'u7^(j»{y;'(j) - 0y2U)} 



n 
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In the case of weak fractional cointegration, we apply Assumption 2.6, the Continuous 
Mapping Theorem and integration by parts to obtain 

n-^£HdA,.,. = -n-'-''i{n)y2weij,n){yUj) - ey^.U)} ^ f h,{t)dZH{t) (4.14) 

where, by independence of B^^ and ^ 

Znit) = c,B^i\\it) - ec,B^S\\,t) ^'=^ Jxl^cl + X^O^cl Bh 



and Bh is a standard fractional Brownian motion. The first part of Lemma 4.4 shows that 
c?Ar/ is negligible under the assumptions of Theorem 3.4. This, and the convergences (4.12), 
(4.13) and (4.14) conclude the proof of Theorem 3.4. □ 

We now prove Theorem 3.5. Since /(.^(O) = h^{l) = 0, we have X]j=o ""^^0' ^) ~ 0' hence 

n—l n—1 n—1 

^wi{j,n)y'^{j) = W£(j,n)(C,,7v,(i) - 6,o) = ^Mj,n)i^i,Ndj) - • 

j=0 j=0 j=0 

Denote Si^ = and for k > I, Si^k = Yl'j=ii^i,N,(j) - f^D- Define uji{j, n) = n{wi{j + l,n)- 
we{j,n)}. Applying again summation by parts, we have 

n—1 -, n—1 



^ we{j, n)y1{i) = -- ^ ^Ki. + Wi{n, n)Si^n-i + Wi{0, n)(^i,o - /^*) 

j=0 j=l 

Under Assumption 2.7, by the Continuous Mapping Theorem, we obtain 

n 

n^~'^£{n)dAyv/ = -n-'^i{n)^Wi{j,n)y^{j) 



4 h';{t)B'£{t) dt - h'{l)B'§{l) ^- f h',{s) di?«(.) . (4.15) 
Jo Jo 

The second part of Lemma 4.4 implies that the term d^r/ does not contribute to the 
limit under the Assumptions of Theorem 3.5. This, and the convergences (4.12), (4.13) 
and (4.15) conclude the proof of Theorem 3.5. □ 

Lemma 4.4. Under the assumptions of Theorem 3.4, then d^,./ = Op{l). Under the 
assumptions of Theorem 3.5, then dAr,e = Op{n~^^'^) . 



Proof. Applying (4.11) to r, we see that we only need to prove that the independence be- 
tween the counting processes and the efficient shocks and (3.3) implies that ^"=1 wi^j 
Op{n) and (3.4) implies that ^"=1 ^^(j, ''^)'"i,j = Op{n^^^). We start with ri. 

n Ni{n) 

i=i k=i ti,fc<i<ti,fc+A2(ti,fc) 
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Taking conditional expectation yields, for g = 1, 2, 



E 



lA/- 



Ni{n) 
k=l 



Applying (3.3) if g = 1 and (3.4) if g = 2 shows that the last term is Op{n). This proves 
that ^"=1 W£(i, 'n.)rij = Op{n) under the assumptions of Theorem 3.4 and Op{^/n) under 
the assumptions of Theorem 3.5. The term ^^j^iWi{jin)r2,j is dealt with similarly. □ 



4.4 Proof of Theorem 3.6 



Write 



with y{i) = yl{j) - Oy'^ij), r{j) = fi(j) - 6'r2(j) and 
ri(s) 

and the DFT is defined as in (3.6). Applying summation by parts as in (4.11), we obtain 



Ni(s) N2{s) 

^ ei,k , r2{s) = ^ e2,k ■ 

fc=A^l(i2,JV2(s)) + l fc=A'2(tl,]Vi(s)) + l 



n-1 



dAy2,£ = y'w^Kj»y2(j) = / Wi{\s/6],n)y2{s)ds = - / Wi{\nt/6],n)y2{ns) dt 

j=0 Jo 

with Wi{j,n) = n{hi{{j + l)/n) — hi{j/n)} as before, and \t~\ is the smallest integer larger 
than or equal to t. This yields 

n-^/^dy,^i 4 -Ee / h'e{s)B{6s) ds = / hi{s) dB{s) . 
Jo Jo 

Since r]j = — ^j-i, we have 

J{j-1)5 

Differencing cancels the term 6^o- Applying (4.11) and summation by parts and the prop- 



erty that X]^=o ""^^0' ''^) = 0, we obtain 



n-l 



- Y] ^iU^ / ^N^is) ds = V Wi{j, n) / {6,7v,(s) - K } ds 

~L J(i-i)s n ~k >^ 0-1)5 



n-l 



n5 



d^xp^. 



j=0 



j=0 



I " ri^ \ i>{n-i)s 

— Y] / {ii,N,{s) - /i- } ds - -w;£(n, n) / {6,JV,(.) - ds 
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Under Assumption 2.8, we thus have, with Z = B^j^^ — 6B^^\ 



n'-^'id^y^^,, - 9d^y.,,} A f K{s)Z{8s) ds - h'{l)Z 

Jo 



(S). 



We must now deal with the remaining terms of the cointegrating error. If if > 1/2, 
Lemma 4.5 imphes that the term (iAf/ does not contribute to the hmit. If if = 1/2, both 
terms are of the same order. This concludes the proof of Theorem 2.8. □ 

Lemma 4.5. Under the assumptions of Theorem 3.6 
Proof. Applying as usual summation by parts, we obtain 

Ni{n5) n „js 



^ k=i ,=1 

k=l 

= ei,fe / wi>{\s/5'\,n)ds 

ei,k / w,{\nt/5'\,n)dt . (4.16) 



k=l J {-ti,k/n)r\S 

Taking conditional expectation and applying (3.4), we obtain 

k=l 



□ 



4.5 Additional Lemmas 

Lemma 4.6. If the durations ti^k — ti^k-i form a stationary ergodic sequence with finite 
moment of order 2p + 1, if P(tj,i > 0) = 1 and if the associated point process has finite 
intensity, then 

supE[(ti,7v^(,)+i - sY] < oo . 

s>0 
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Proof of Lemma 4-6. We omit the index i. Let 6t denote the shift operator and let A{t) 
be the forward recurrence time. Then A{s) = tAr(s)+i — s = ti o Og. Since the sequence 
{tj} is stationary under P, there exists a probabihty law P* such that is a stationary 
ergodic point process under P*, see BacccUi and Brcmaud (2003, Section 1.3.5). Applying 
Baccelli and Bremaud (2003, Formula 1.3.3), we obtain 



E[AP{s)] = X-^E* 



N{1) 

k=l 
■jv(i) 

k=l 



i-s- 



N{1) 



k=l 



< X~^E* 



N{1) 



sY 



k=l 



A-^r [iV(l){tjv(s+i)+i - sY] < A-i{E*[Ar(l)2]}V2|E*[(t 



Af(s+1)+1 



2pj|l/2 



(4.17) 



Since is stationary under P*, the last term does not depend on s, and by the Ryll- 
Nardzewski inversion formula (Baccelli and Bremaud (2003, Formula 1.2.25)), we have 

r[(tAr(,+i)+i - sfP] = E*[{h + l)'n = AE[ f\h + 1 - s)^Pds < XE[{l+h)^P+'] 

Jo 

By Baccelli and Brcmaud (2003, Property 1.6.3), the point process A^ is stationary and 
ergodic under P* since the sequence of durations is stationary and ergodic. Thus, By 
Daley and Vere-Jones (2003, Theorem 3. 5. Ill), E*[A^(0, 1)^] < oo. Plugging the last two 
bounds into (4.17), we obtain that E[A^(s)] is uniformly bounded. □ 



Lemma 4.7. Assume that there exists an increasing sequence {sn,n > 0} such that sq = 
and 



(a) f is either constant or strictly increasing and differentiable on (s„, s„_|_i) and the jumps 
of f occur at some (but not necessarily all) of the 

(b) if f is eiter constant or increasing on both intervals (s„, Sn+i) and {Sn+i, Sn+2), then f 
has a jump at Sn+i- 



Assume moreover that 



• (minimum duration of trading and nontrading periods) there exists 60 > such that 
Sn+i — Sn> So for all n > 0; 

• (maximum duration of nontrading periods) there exists Cq such that for all n > 0, if 

f is constant on (s„, Sn+i), then Sn+i — s„ < Cq; 

• (non stoppage of time during trading periods) there exists 61 > such that for all 
n >0, f is either constant on {sn, Sn+i), or f'{t) > 61 for all t G (s^, Sn+ij- 
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Let N be a point process with event times {tk} and let N be the point process de- 
fined by N{-) = N{f{-)) with event times {tk}- // supg>o lE[(t;v(5)4.]^ — sY] < oo, then 

SUp^>olE[(tAr(s)+i - sY] < OO. 

Proof of Lemma 4- 7. Define the nondecreasing left-continuous inverse of a nondecreas- 
ing cadlag function / by 

f^iu)=mi{t\fit)>u}. 

Note first that /^(u) < t if only if n < f{t) and /^(/(t)) < t. Thus we see that 

f^{in) <t^L< fit) 

^ Nifit)) > n 
^ N{t) > n . 

This characterizes the sequence {tn}, thus we obtain that t„ = f^{tn). The assumptions 
on / imply the following properties of 

• The jumps of correspond to the intervals (s„, s„+i) where / is constant. More 
precisely, if / is constant on (s„, Sn+i), then has a jump at of size — s„. 
Since is left continuous, it holds that 

f^ifiSn)) = Sn , lim = Sn+1 ■ 

Thus the jumps of are of size Cq at most. 

• If / is increasing on an interval (s^, then is differentiable on f{s~)) 
and < for all t G (/(s^, /(s"))- 

• The jumps of / create no singularity in f*^. If > f{s~), then is constant 
on the interval (/(s"), 

Let [x] denote the smallest integer greater than or equal to the real number x. Then, for 

< s < t, 

< f^{t) - < Co + S^\t - s) . 

Thus, there exits constants ci, C2 such that for all s < t, 

0<f{t)-f{s)<c, + C2{t-s). 
Consider now the forward recurrence time of the point process A^. Then 
< - s = r - r (/(s)) + r (/(.)) - s 
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Thus, there exists constants C3 and C4 such that 

supE[(t7v(s)+i — sY] < C3 + C4SupE[(t 



s>0 



s>0 



N{s)+1 



□ 

Lemma 4.8. Let {et} be a sequence of i.i.d. positive random variables with finite mean 
H^. Let {Yk} be a stationary standard Gaussian process such that 



for H G (1/2, 1) and I a slowly varying function. For k > 1, define 



(4.18) 



Then the sequence {rk} is ergodic and Assumption 2.1 holds with \ = fi^e'^^^^. //P(ei > 
0) = 1 the Assumption 2.2 holds with jj, = \ = iJi~^e~" // moreover ]E[e^] < 00 for all 
q > 1, then (3.3) and (3.4) hold. 

Remark 4.1. If instead of (4.18) we assume that 

00 

|cov(Fo,'^fc)| < 00 , 

k=l 

then the moment requirement can be relaxed to IE[ei] < 00 to obtain (3.3) and IE[ef] < 00 
to obtain (3.4). 

Proof of Lemma 4-8. Note first that IE[r^] < 00 as long as IE[ei] < 00. By Lemma 4.6, in 
order to check condition (3.3), we must only prove that the induced point process has finite 
intensiy, i.e. there exists t > such that E[iV(t)] < 00. See Baccelli and Bremaud (2003, 
Section 1.3.5). Note that 

00 00 
E[N{x)] =^F{N{x) >k) = ^P(tfc < x) . 



k=l 



k=l 



Thus, it suffices to prove that the series on the righthand side is summable. Denote 
fj, = E[rfc] and p„ = cov(Yo,Yn). Applying Deo et al. (2009, Proposition 1), we have 



E 



^ Tfc - n/i 

k=l 



PI 



with Vn = n^i{n). If E[e^] < 00 for p such that p{l — H) > 1, for n such that nfj, > x, it 
holds that 

ntk <x) = oix-'vi) 



and this series is summable. 



□ 
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Lemma 4.9. Assume that {Tk} and {C,k} are mutually independent stationary sequences 
such that E[^fc] = 0, lE[r|] < oo and IE[^^] < oo. Assume that the sequence of durations is 
weakly stationary and that cov(ro, r^) = 0(n~'') for some S > and supg>o IE[tAr(s)+i — s] < 
oo. Assume that cov(,^i,,^„) ~ cn"^^'^, with H G (1/2, 1) and c> 0, and that 

k=l 

for some c' > 0. Then 

Jo 

for some c" > 0. 



Proof of Lemma 4-9. Denote E[rfc] = /i > 0. 



/ ^N{s) ds — ^ Tk+lC,k - {tN{T)+l - T)C,iy{T)+l 
JO 1,-n 



k=0 

N{T) N(T) 



k=0 k=0 



By independence of {rk} and {C,k}, we have (assuming without loss of generahty that 
2H-6> 1), 



var ^(r^+i - /i)^fc = 0{n 



2H~S\ 



.fc=0 



Thus, Ylt^^oi'^k+i — /^)^fc =^ 0. Hence by the continuous mapping theorem, it also 
holds that ^k^ \'''k+i ~ /^)'Cfc ^ 0. By independence and by assumption, {^fq^t-^j^x — 
T)^N(T) = Op{l). By the continuous mapping theorem, J2k^^^ =^ d Bh{^~H). □ 

Lemma 4.10. Let {rk}, {Vk} and {C,k} he sequences of random variables such that 



• {Cfc} an i.i.d. sequence of zero-mean and unit variance random variables; {r^} and 
{Vfc} are sequences of positive random variables; 

• the sequences {(t^, V^)} and {C^k} are mutually independent; 

• there exists s > such that X]fc=i '^k+i^k ~^ ^'^i 

• supfc>oE[r^^+^V^^+''] < oo for some e > 0; 

• sup^>oE[t7v(s)+i - s] < oo. 
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Define = CkVk- Then T ^1'^ ^Af(s) ds ^ cB for some c > 0. 

Proof. Let J-^ be the sigma- field generated by random variables {tj+i, (j, Vj,j < k}. Then 
E[,^fcrfc_|_i I J^k-i] = Tfc+iVfcE[Cfc] = 0. Thus, 1x^+1,^^} is a martingale difference sequence. 
Under the stated assumptions, the martingale invariance principle Hall and Hcydc (1980, 
Theorem 4.1) yields that n~^^^ Yl^k=i'^k+i^k =^ cB for some c > 0. As in the proof of 
Lemma 4.9, denote ]E[rfc] = /i > and write 

/ iN{s) ds = ^ Tkj^lik + (^iV(T)+l - T)^M(T) . 
•^0 fc=0 

By the continuous mapping theorem, we have that T~^/^ Xlfc^T ^ '^k^-i =^ Aci?. As previ- 
ously, the last term is a neghgible edge effect. This concludes the proof. □ 
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